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We present a search for the lightest supersymmetric partner of the top quark in proton-antiproton 
collisions at a center-of-mass energy >/i = 1.96 TeV. This search was conducted within the framework 
of the imparity conserving minimal supersymmetric extension of the standard model, assuming the 
stop decays dominant ly to a lepton, a sneutrino, and a bottom quark. We searched for events 
with two oppositely-charged leptons, at least one jet, and missing transverse energy in a data 
sample corresponding to an integrated luminosity of 1 fb _1 collected by the CDF experiment. No 
significant evidence of a stop quark signal was found. Exclusion limits at 95% confidence level in the 
stop quark versus sneutrino mass plane are set. Stop quark masses up to 180 GeV/c 2 are excluded 
for sneutrino masses around 45 GeV/c 2 , and sneutrino masses up to 116 GeV/c 2 are excluded for 
stop quark masses around 150 GeV/c 2 . 
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The minimal supersymmetric standard model [1] 
(MSSM) was introduced to solve several problems that 
arise in the standard model (SM). These include: the 
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hierarchy problem that requires the fine-tuning of theo- 
retical parameters in order to obtain cancellation of large 
quantum corrections to the Higgs mass; the lack of con- 
vergence of the strong, weak, and electromagnetic gauge 
couplings at the grand-unification energy scale; and the 
lack of a dark matter candidate. 

The MSSM assigns a new bosonic counterpart to each 
SM fermion and likewise a fermionic superpartner to each 
SM boson. This results in scalar partners and qR to 
the SM quark helicity states qL and qR. There can be 
two supersymmetric mass eigenstates for each supersym- 
metric quark (squark) corresponding to the two fermionic 
degrees-of-freedom of the SM quark. The supersymmet- 
ric scalar top quark (stop) mass eigenstates t\ and £2 
are rotated relative to II and Ir by a mixing angle 0$. 
In some models [2|, 2 ~ m t . The large mass of the 
top quark and the corresponding large value of the top- 
to-Higgs coupling constant may lead to a large splitting 
between and m^ 2 . Consequently the lower mass stop 
quark eigenstate is expected to be the lightest of all the 
squarks, with a mass even below the top quark, making 
its detection at the Tevatron a realistic possibility. 

The MSSM possesses a new conserved quantity called 
i?-parity (i? p ), defined as R p = (_i)3(b-l)+2S ? where 
B is the baryon number, L is the lepton number, and S 
is the spin. As a consequence, the lightest supersymmet- 
ric particle (LSP) must be stable, and is a dark matter 
candidate. Because the initial state of pp collisions has 
R p = +1 and supersymmetric particles have R p = — 1, 
supersymmetric particles must be pair-produced. At 
the Tevatron stop quarks are expected to be produced 
primarily through gluon-gluon (gg) fusion and quark- 
antiquark (qq) annihilation, with gg fusion dominant at 
low stop masses (< 100 GeV/cf) and qq annihilation 
dominant at higher stop masses [3j. 

The produced stop (t) quarks can decay via several 
possible channels, depending on the masses of the par- 
ticles involved. Two-body decays include t — >> txi, 
t bxi , and t — >> cxi where Xi 1S the lightest neu- 
tralino and xt the lightest chargino. These decays may 
not be kinematically possible for a light stop or they may 
be suppressed by higher order diagrams. They are also 
constrained by existing limits 0, Hj . Possible three-body 
decays include t — >> W+bxi, t — » bl+vi, and t — >> bl+vi 
where I is the supersymmetric lepton and v\ the super- 
symmetric neutrino. Limits on the supersymmetric lep- 
ton and neutralino masses from experiments at the large 
electron-positron collider (LEP) [4[ restrict the range of 
stop masses available to the first two decay modes. 

The decay t — >> bl + i>i, which proceeds via a virtual 
chargino, is the subject of this analysis. We assume the 
branching ratio for this decay mode is 100% and that 
electrons, muons, and taus are equally likely decay prod- 
ucts. While electrons and muons are detected directly, 
taus are only included opportunistically in this analysis 
through their decays into electrons and muons. We also 
assume that the supersymmetric neutrino decays neu- 
trally into the LSP (or is the LSP), thus escaping unde- 



tected and leading to missing transverse energy (J£t) 
in the detector. Since stop quarks are produced in pairs, 
we search for events with two opposite-charge leptons (ee, 
e/i, /i/i), ^t, and at least one hadronic jet. 

Previous searches at LEP and at the Tevatron 0, Q 
for the stop quark using the same topology as this analy- 
sis have produced negative results. These analyses have 
set 95% confidence level exclusion limits in the stop- 
sneutrino mass plane. This paper extends the earlier 
CDF results based on 107 pb" 1 of data at = 1.8 TeV 
to 1 fb- 1 at v^= 1-96 TeV. 

The structure of this paper is as follows: Section II de- 
tails the detector and the data set; Sections III describes 
the background estimation; Section IV the signal pre- 
dictions and systematic uncertainties; Section V event 
pre-selection and control samples; Section VI explains 
the optimization of event selection cuts; and Section VII 
presents the results and conclusions. 

II. THE DETECTOR AND DATA SET 

The data sample corresponds to an integrated luminos- 
ity of 1 fb -1 of pp collisions at a center-of-mass energy of 
1.96 TeV collected with the CDF detector [7] at the Fer- 
milab Tevatron. Of particular relevance to this analysis 
are the tracking system, the calorimetry, and the muon 
detectors. 

The tracking system consists of two silicon micro-strip 
detectors and a multi-wire open-cell drift chamber, the 
central outer tracker (COT). The silicon vertex detec- 
tor (SVX) and the intermediate silicon layers cover the 
pseudorapidity [6[ region \rj\ < 2 while the COT covers 
1 77 1 < 1. The tracking system is surrounded by a su- 
perconducting solenoid with a magnetic field of 1.4 T. 
The relative track momentum resolution provided by the 
COT is a PT /p T 2 ~ 0.0015(GeV/c)- 1 . 

Outside the magnet in the radial direction are electro- 
magnetic and hadronic calorimeters arranged in a pro- 
jective tower geometry with a tower granularity of At] x 
Acj) c± 0.1 x 0.26 in the central region. The central elec- 
tromagnetic calorimeter (CEM) utilizes lead-scintillator 
sampling and covers \rj\ < 1.1, with energy resolution 
for electrons (j Et /E t ~ 13.5%/ (GeV) © 2%. The 
central electromagnetic shower maximum detector (CES) 
is located inside the CEM at a depth of six radiation 
lengths, close to the position of maximum electromag- 
netic shower development, and is used for the determi- 
nation of the shower shape and for matching the shower 
location to the track extrapolation. The central hadron 
calorimeter (CHA) uses iron-scintillator sampling and 
covers \r]\ < 0.9. Its energy resolution for hadrons is 
a ET /E T ~ 75%/V^t (GeV) ©3%. 

Additional calorimetry extends the coverage in the 
forward direction and is used in this analysis for cal- 
culating and jet energies but not for lepton iden- 
tification. The plug electromagnetic calorimeter covers 
1.1 < 1 77 1 < 3.6 and is constructed of lead and scintil- 
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FIG. 1: Pt distributions in the pre-signal region for (a) the ee channel, (b) the e\i channel, (c) the channel, and (d) the 
three channels combined. Data are shown as the points with error bars (statistical only). Shown as stacked histograms are 
the backgrounds arising from misidentified hadrons and decays-in-flight (7+fake), bb and cc (HF), DY, dibosons, and tt. For 
reference, the expected signal for (m^rriu) = (150,75) GeV/c 2 , multiplied by five, is shown as the dashed line. 



lator layers with an energy resolution for electrons of 
ct Et /E t ~ 16%/V^r (GeV) © 1%. The plug hadron 
calorimeter covers 1.3 < \rj\ < 3.6 and is constructed 
of iron and scintillator layers wit h an energy resolution 
for hadrons of ct Et /E t ~ 80%/V^r (GeV) © 5%. The 
iron and scintillator wall hadron calorimeter covers the 
intermediate region 0.7 < |^| < 1.3 between the central 
and plug hadron calorimeters. Its energy resolution is 
ct Et /E t ~ 75%/V^t (GeV) © 4%. 

Radially outside the calorimetry is the muon detection 
system. The parts of the muon detector relevant to this 
analysis are the central muon detector (CMU), the cen- 
tral muon upgrade detector (CMP) and the central muon 
extension (CMX). The CMU consists of four-layer drift 
chambers and covers the pseudorapidity range |^| < 0.6. 
The CMP is made of four layers of single-wire drift cells 
located behind an additional 3.3 interaction lengths of 
steel and covers \rj\ < 0.6. Muons with reconstructed 
track stubs found in both the CMU and CMP are la- 
beled CMUP muons. The CMX extends the coverage to 
1 77 1 < 1 and is made up of eight layers of drift tubes. 



The data were collected using a three-level trigger sys- 
tem. The first two levels are mostly hardware-based 
while the third level is software-based and is a fast ver- 
sion of the offline event reconstruction package. The 
online selection requires at least two lepton candidates 
falling into the following categories: CEM-CEM, CEM- 
CMUP, CEM-CMX, CMUP-CMUP, and CMUP-CMX, 
where the leptons are labeled by the detector compo- 
nents used in their identification. The triggers used in 
this analysis had an Et threshold of 4 GeV for electrons 
and a pt threshold of 4 GeV j c for muons. 

All lepton candidates require the presence of a well- 
reconstructed track in the COT detector. Offline, elec- 
trons are required to pass a x 2 comparison of the CES 
lateral shower profile in the r — z view and the profile 
extracted from electron test-beam data. The shower po- 
sition in the CES must also match the extrapolated po- 
sition of the track. In addition, the lateral shower pro- 
file in the CEM towers must be consistent with that ex- 
pected from test-beam data. The energy deposited in 
the calorimeters must be consistent with the track mo- 
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FIG. 2: Ht distributions in the pre-signal region, shown as in Fig. \T\ 



mentum measurement and the energy deposition in the 
CHA must be small. To reduce the background arising 
from the decay of hadrons in jets, the electrons are re- 
quired to be isolated. Isolated electrons are selected by 
requiring that the remaining transverse energy after sub- 
tracting the transverse energy associated with the elec- 
tron in a cone of 0.4 in rj — (j) space (E™°) must be less 
than 2 GeV if E T < 20 GeV or E^°/E T be less than 
0.1 if Et >20 GeV. Electrons arising from photons that 
convert into e + e~ pairs are removed by cuts applied to 
pairs of opposite-charge tracks with small opening an- 
gles. The electron identification efficiency was measured 
with Drell-Yan (DY) electrons to range from 75 to 83% 
[8] increasing with the electron E^. 

Muons are identified by matching tracks with recon- 
structed track stubs in the CMUP or CMX. The energy 
deposited in the calorimeters must be consistent with a 
minimum- ionizing particle. The isolation requirement is 
similar to that for electrons: E™° < 2 GeV if pt < 20 
GeV/c or / 'E T < 0.1 if p T > 20 GeV/c. The muon 
identification efficiency was measured with J/ip and Z 
data to be in the range 90 — 96% [H rising with the muon 
PT- 

For both electrons and muons the impact parameter 



(do) of the track with respect to the beam line position 
must be less than 0.2 cm if the track is based on COT 
information only or less than 0.02 cm if silicon-based 
tracking is also available. The longitudinal position of 
the event vertex is required to be within 60 cm of the 
center of the detector. 

Jets are identified in r] — <p space as a group of elec- 
tromagnetic and hadronic towers using a clustering algo- 
rithm with cone size AR = \/(Ar]) 2 + (A</>) 2 = 0.7. The 
jet energy resolution is a ^ 0.1 • Et (GeV) ® 1 GeV for 
jets with 35 < E T < 450 GeV. In order to find genuine 
hadronic jets and to avoid counting electrons and pho- 
tons as jets, the fraction of the total energy deposited 
in the electromagnetic calorimeters is required to be be- 
tween 0.1 and 0.9. Corrections to the jet energy are ap- 
plied to take into account 77-dependent losses, luminosity- 
dependent multiple interactions, and non-linearities in 
calorimeter response @. In this analysis we require at 
least one jet with \rj\ < 2 and E T > 15 GeV after correc- 
tions. 

is calculated as the magnitude of the negative 
vector sum of the transverse energy deposited in the 
calorimeter towers with |^| < 3.6 and energy larger than 
0.1 GeV. Since muons are minimum- ionizing particles 
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FIG. 3: pt distributions of the highest pt lepton in the pre-signal region, shown as in Fig. [T] 



and deposit only ~ 2 GeV energy in the calorimetry, 
the muon momentum from tracking (minus the expected 
energy deposition) is used to correct fJ^, The energy 
corrections applied to jets with \rj\ < 2.4 and E^ > 10 
GeV are also propagated to the calculation of ^t- 



III. BACKGROUND ESTIMATION 

Several background sources result in events with dilep- 
tons, jets, and These backgrounds are tt production, 
other heavy-flavor quark (66 and cc) production, DY pro- 
duction of lepton pairs where mismeasurement of a jet 
or lepton results in substantial ^t, diboson production 
(WW, WZ, ZZ, and W7), and events with a lepton 
and a misidentified or secondary lepton (Z+fake). The 
jets in most backgrounds result from QCD initial-state 
or final-state radiation (ISR or FSR). 

The tt, DY, and diboson backgrounds are estimated by 
generating Monte Carlo (MC) events using the pythia 
[10^ event generator followed by a GEANT-based pj| de- 
tector simulation. Backgrounds arising from fake leptons 
and heavy flavor (66, cc) are estimated using data-driven 
methods. 



The MC tt samples are normalized to the next-to- 
leading-order (NLO) cross section a t i = 6.71 pb [12] with 
a top mass of 175 GeV/c 2 . The Z — >• ee, Z jiji, and 
Z — >• tt samples are normalized to a leading-order cross 
section times branching ratio of 1272 pb for m\\ > 10 
GeV/c 2 times a k- factor of 1.4 [13[ to correct for NLO 
contributions. The NLO cross sections times branching 
ratios for WW llvv, WZ -> IIX, and ZZ 11 X 
are 1.27, 0.365, and 1.513 pb respectively For 
W7 —> evy and W7 —> [ivy, we use the leading order 
cross section of 21.5 pb times a fc- factor of 1.34 to ac- 
count for NLO contributions [l5| ]. 

The lepton plus fake background consists of events 
with a genuine lepton plus a "fake" lepton which is either 
a light hadron misidentified as a lepton or an uninterest- 
ing lepton from pion or kaon decay-in-flight. For muons 
the fakes can be particles that penetrate the calorimeters 
and absorbing material and reach the muon detectors or 
decay in flight to muons. In the case of electrons, fakes 
are usually jets that are misidentified as electrons, mainly 
due to neutral pions that decay to photons which shower 
in the electromagnetic calorimeters. This background is 
estimated by examining samples of single-lepton events 
taken with the Et > 8 GeV electron calibration trigger 
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FIG. 4: pt distributions of the next-to-highest pt lepton in the pre-signal region, shown as in Fig. \T\ 



and the pt > 8 GeV/c muon calibration trigger. Events 
with at least one E T > 4 GeV central jet (\rj\ < 1.1) or 
one isolated track with p^ > 4 GeV/c passing the muon 
track quality cuts, excluding the trigger lepton, form the 
/+fake-electron and /+fake-muon candidate samples, re- 
spectively. To determine the background, it is necessary 
to estimate a fake rate from samples of jet events that 
contain a negligible number of directly produced leptons. 
Four samples containing jets with Et > 20, 50, 70, and 
100 GeV, respectively, are used. The electron fake rate 
is defined as the probability of a jet being misidentified 
as an electron and the muon fake rate as the probabil- 
ity of an isolated track being misidentified as a muon. 
The fake rates are determined as a function of jet Et 
for electrons and of track px for muons. The fake rate 
for 20 GeV electrons is about 0.0002 fakes per jet and 
the fake rate for 20 GeV/c muons is about 0.004 fakes 
per track. These fake rates are then applied to each 
fake-electron and fake- muon candidate, one at a time, in 
the single lepton sample. Events that pass the analysis 
cuts are assigned a weight, which is the appropriate fake 
rate, scaled to the integrated luminosity of the dilepton 
data sample relative to the integrated luminosity of the 
single lepton sample, taking into account the single lep- 



ton trigger prescale. The sum of the weights then forms 
the /+fake background. A 50% systematic uncertainty 
on the misidentifed lepton background is assigned based 
upon the differences between the four jet samples used 
to determine the fake rates. 

The background arising from heavy-flavor (66, cc) is 
estimated utilizing dilepton events enriched in heavy- 
flavor events by inverting the normal impact parameter 
requirements, i.e. \do\ > 0.2 cm for COT-only tracks and 
| do | > 0.02 cm if the track included information from 
the SVX, and requiring that at least one lepton pass the 
inverse do cuts. A "scaling region" is defined as the dilep- 
ton invariant mass range 15 < M\\ < 35 GeV/c 2 . The 
scaling region has no requirements imposed on the 
number of jets, or other kinematic variables. In this re- 
gion the only significant contributions to the data sam- 
ple with normal do cuts are due to DY, heavy-flavor, and 
light hadrons misidentified as leptons. The DY contri- 
bution is derived from MC samples, and the misidenti- 
fied lepton component is estimated using the technique 
described above. The remaining events are attributed 
to heavy-flavor (HF) and are used to calculate scaling 
factors, defined as the ratio of these inferred HF events 
passing the normal do cuts to the total number of events 
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FIG. 5: The cut group b pT distributions for (a) the ee channel, (b) the e\i channel, (c) the fifi channel, and (d) the three 
channels combined with all other cuts applied. Data are shown as the points with error bars (statistical only). Shown as stacked 
histograms are the backgrounds arising from misidentified hadrons and decays-in-flight (7+fake), bb and cc (HF), DY, dibosons, 
and it. MC predictions for (m^rriu) = (140,90) GeV/c 2 are stacked on top of the sums of backgrounds and are shown as the 
dashed line. The cut selects events to the right of the vertical line. 



passing the above "inverse do" cuts. It is assumed that 
this scaling ratio is independent of dilepton mass in re- 
gions not dominated by Z boson decays. To estimate 
the heavy-flavor background, the resulting scaling fac- 
tors, typically of order 2, are then applied to "inverse 
do" events passing cuts appropriate to the various con- 
trol regions as well as the signal region. No heavy-flavor 
contribution to the signal region survives our final cuts. 
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IV. SIGNAL SAMPLES AND SYSTEMATIC 
UNCERTAINTIES 

A total of 74 MC signal samples corresponding to 
values of m~ t = [55,190] GeV/c 2 and = [45,110] 
GeV/c 2 were generated using pythia and the geant- 
based detector simulation. The events are scaled using 
the NLO PROSPINO cross section [H, [13 calculated with 
the CTEQ6M [18] parton distribution functions (PDFs). 
The cross sections depend on the stop quark mass and 



TABLE I: Expected backgrounds (DY, tt, diboson, misiden- 
tified hadrons and decays-in-flight (Z+fake), and bb and cc 
(HF)) and observed events in the pre-signal region. The 
first uncertainty listed is statistical and the second system- 
atic. The systematic uncertainty includes a 6% uncertainty 
on the luminosity common to all entries. 
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FIG. 6: The cut group b distributions of the absolute value of the azimuthal angle between the dilepton system and ^t, shown 
as in Fig. [5] 
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TABLE II: Expected DY background and observed events 
in the Z control region. The first uncertainty listed is statis- 
tical and the second systematic. The systematic uncertainty 
includes a 6% uncertainty on the luminosity common to both 
channels. 
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only weakly, through higher-order corrections, on other 
SUSY parameters. Limits on the production cross sec- 
tion can therefore be translated into lower limits on the 
lightest stop quark mass without reference to other SUSY 
parameters. The branching ratio for the decay t — >• 
where I = e, \i or r with equal probablity, is assumed to 
be 100%. 

There are several sources of systematic uncertainty on 
the background estimation. The systematic uncertainty 
due to the jet energy scale is determined by varying the 
jet energy corrections by ±la |9]. The resulting uncer- 
tainty varies from less than 1% to 35% and is largest 



TABLE III: Expected backgrounds (DY, ti, diboson, 
misidentifled hadrons and decays-in-flight (e+fake, /i+fake), 
and bb and cc (HF)) and observed events for the high-^x / 
no jet control region. The first uncertainty listed is statisti- 
cal and the second systematic. The systematic uncertainty 
includes a 6% uncertainty on the luminosity common to all 
entries. 
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FIG. 7: The cut group b pT2 distributions, shown as in Fig. [5] 
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TABLE IV: Expected backgrounds (DY, tt, diboson, 
misidentified hadrons and decays-in-flight (e+fake, /i+fake), 
and bb and cc (HF)) and observed events for the 1ow-^t / 
no jet control region. The first uncertainty listed is statisti- 
cal and the second systematic. The systematic uncertainty 
includes a 6% uncertainty on the luminosity common to all 
entries. 



TABLE V: Expected backgrounds (DY, tf, diboson, misiden- 
tified hadrons and decays-in-flight (e+fake, /i+fake), and bb 
and cc (HF)) and observed events for the \ow-pT / one or 
more jet control region. The first uncertainty listed is statis- 
tical and the second systematic. The systematic uncertainty 
includes a 6% uncertainty on the luminosity common to all 
entries. 
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FIG. 8: The cut group b Ht distributions, shown as in Fig. [5] The cut selects events between the vertical lines. 
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TABLE VI: Table of cuts for the cut groups defining the four 
signal regions. 



ficiency, 2% on lepton identification efficiency, and 50% 
on the number of misidentified electrons and muons. The 
uncertainties on the cross sections used in the MC gen- 
eration of the background are 8% for ti, 2% for DY, 6% 
for WW, 8% for WZ, 10% for ZZ, and 7% for Wj. 

The systematic uncertainties on the MC stop signal 
estimation are nearly identical to those on the back- 
ground estimation. The uncertainty arising from the jet 
energy scale varies from 1% to 11% depending on the 
stop-sneutrino mass difference. The uncertainties due to 
ISR/FSR, PDFs, luminosity, trigger efficiency, and lep- 
ton identification are the same as for the background es- 
timation. 



for DY events, which typically contain low-i^x jets. The 
uncertainty arising from ISR and FSR is determined by 
varying the parameters in PYTHIA that control the gen- 
eration of ISR/FSR. The resulting uncertainty is c± 3%. 
The uncertainty on the acceptance arising from the PDFs 
used in the MC is estimated using the uncertainties on 
the CTEQ eigenvectors and determined to be 2%. Other 
systematic uncertainties are: 6% on the measurement of 
the integrated luminosity, 2% on the dilepton trigger ef- 



V. INITIAL EVENT SELECTION AND 
CONTROL SAMPLES 

We first define a pre-signal region by applying several 
event selection cuts to significantly reduce background 
and provide a data sample loosely consistent with the 
stop quark signature. We subsequently optimize addi- 
tional cuts to improve the sensitivity of the search. This 
is done prior to revealing the contents of the data in 



13 




20 40 60 80 100 120 2 140 

m T (GeV/c ) 




20 40 60 80 100 120 2 140 

m T (GeV/c ) 




20 40 60 80 100 120 2 140 

m T (GeV/c ) 




20 40 60 80 100 120 2 140 

m T (GeV/c ) 



FIG. 9: The cut group b mi distributions, shown as in Fig. [5] 



the pre-signal region. At the pre-signal stage the follow- 
ing cut requirements are applied: two opposite-charge 
leptons, one with Pt> 10 GeV/c and the other with 
Pt> 5 GeV/c; mu > 15 GeV/c 2 , in order to remove 
sequential 5-hadron decays and low mass resonances; 
mu < 76 GeV/c 2 or mu > 106 GeV/c 2 for same-flavor 
dilepton events, in order to eliminate Z boson events; at 
least one jet with corrected Et > 15 GeV and \rj\ < 2; 
^ T > 15 GeV; A#(e,highest-£;T jet) > 0.4; AR(l,l) > 
0.4; and > 20° between JS^ and each of the leptons 
and the hightest Et jet. 

Figure [T] shows the ]£t distributions separately for ee, 
e/i and fifi events, as well as the summed distribution, 
in the pre-signal region. The expected $t distribution 
for stop quark events with stop mass 150 GeV/c 2 and 
sneutrino mass 75 GeV/c 2 is also shown scaled up by a 
factor of five. Figures [2]H] show the corresponding plots 
for i?T 5 the pt of the highest pt lepton (pti) 5 and the pt 
of the next-to-highest pt lepton (pt2)- Ht is defined as 
Ht = -^T+PTi+PT2+^Tj, where ^Tj is the transverse 
energy of the highest pt jet. Table U lists the sources 
of expected background for the pre-signal region and the 
number of observed events. Good agreement of data with 
the background estimations is observed. 



To check the accuracy of our estimation of SM back- 
grounds, a number of "control regions" are denned. One 
control region consists of same-flavor, opposite-charge 
lepton events with invariant mass 76 < mu < 106 GeV/c 2 
(the Z region). No jet or requirements are imposed. 
The background is expected to arise almost entirely from 
DY processes. Table [XT] gives the number of observed 
events in the Z region and the expected background. 
Good agreement is seen, demonstrating accurate mod- 
eling of the predominantly DY background. 

Another control region consists of events where we re- 
quire ]S T > 15 GeV and no jets with E T > 15 GeV. All the 
other pre-signal cuts are applied except for the Ac/) > 20° 
cut between and each of the two highest pt leptons 
and the highest E^ jet. This control region predom- 
inantly tests the modeling of the high tail of DY 
and HF events. Table EH gives a breakdown of the ex- 
pected backgrounds and the number of observed events 
for this control region. Reasonable agreement is observed 
between the expected background and data. The signifi- 
cant number of e\i DY events (comparable to the number 
of ee and fifi events) is due to Z —> rr decays where one 
r subsequently decays to an electron and the other to a 
muon. 
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95% CL upper limit on stop pair production 
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FIG. 10: Stop cross section upper limits for fixed stop-sneutrino mass differences. The NLO prospino cross section with 
CTEQ6M PDFs is shown as the solid curve. The band represents the theoretical uncertainty in the cross section due to 
uncertainties on the renormalization and factorization scales and the PDFs. 



Two additional control regions are defined as having 
< 15 GeV and either the presence or absence of jets 
with Et > 15 GeV. These control regions are sensitive 
to the modeling of heavy flavor and the photon compo- 
nent of DY. All other cuts are the same as for the pre- 
vious control region. Tables [IV] and [V] list the expected 
backgrounds and the number of observed events for these 
control regions. In both cases, the data agree well with 
the expected backgrounds. 

VI. OPTIMIZATION OF CUTS 

The signal regions are determined by optimizing cuts 
on variables that minimize the limit on the it pro- 
duction cross section determined with Bayesian meth- 



ods [l9[. The ee,e/i, and fifi final states are tuned si- 
multaneously and flavor-independent cut values are ob- 
tained. The most useful variables for discriminating stop 
quark signal from background are: J£t] Ht, the abso- 
lute value of the difference in azimuthal angle between 
the dilepton system and ^t; Pt2] the transverse mo- 
mentum of the dilepton system (p 1 ^); and the trans- 
verse mass between each lepton (1) and ^t, defined as 
mxi = v2pTi^r[l — cos(A(/>(^, j£r))] where the index i 
labels the lepton. 

The optimum values for the cut variables depend on 
the mass difference (Am = m^ — rriy) which we group into 
four sets labeled a through d in Am bands parallel to and 
below the m^ = m^+m^ kinematic limit. Four cut groups 
were chosen as a compromise between having the cuts 
near their optimal values for each stop-sneutrino point 



Scalar top in the dilepton channel: t -» I v b 
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FIG. 11: Observed and expected limits in the stop-sneutrino mass plane. The LEP limits are for a mixing angle of zero which 
provides the greatest reach. 



while keeping the number of cut groups to a minimum. 
The definitions of the four cut groups and the values of 
the cuts used are given in Table fVTl 

The distributions for the signal region are shown 
in Fig. [5] for cut group 6, where the other cuts listed in 
Table [VTl are all applied. The individual backgrounds are 
shown as well as the data. For reference the expected sig- 
nal from the stop-sneutrino mass point (140,90) GeV/c 2 
is also shown added to the stacked backgrounds. The ver- 
tical line represents the lower bound placed on for this 
cut group. The ft? cut is the most effective at reducing 
the DY, bb and cc, and /+fake backgrounds. In particu- 
lar, the DY background is expected to show a small en- 
ergy imbalance in the detector, while the reference-point 
signal events are characterized by large energy imbalance. 
The cut on fl^ reduces the DY background by more than 
a factor of 10 even when all other final cuts have been 



applied, while reducing the expected signal by only 20%. 

As shown in Fig. [6j the cut on the absolute value of 
the azimuthal angle between the dilepton system and J£t 
is useful in further supressing the dominant DY back- 
ground. This cut discriminates against DY events where 
the arises from mismeasurement of the leptons or 
from tt events where neutrinos from the r decay result 
in aligned with the leptons coming from the decay. 

Leptons arising from stop quark decay are typically 
more energetic than those coming from b or c quark de- 
cay because of the higher stop quark mass. A cut on 
PT2, whose distribution is shown in Fig. [7J removes all 
the bb and cc background remaining after the other cuts 
have been applied. Only 5% of the reference-point signal 
events are removed by this cut. 

The tt background is especially difficult to reduce with- 
out severely impacting the efficiency for stop quark detec- 
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TABLE VII: Expected backgrounds and observed events for cut groups a - d. The uncertainties represent the statistical and 
systematic uncertainties added in quadrature. The systematic uncertainty includes a 6% uncertainty on the luminosity common 
to all entries. Also shown are the number of expected tt events for (stop, sneutrino) masses (130,95), (140,90), (150,75), and 
(180,50) GeV/c 2 in cut groups a, b, c and d respectively. 



tion. However, because the sneutrino carries a significant 
fraction of the available energy, the remaining stop quark 
decay products typically have less energy than their coun- 
terparts from top quark decay. As a result the H T distri- 
bution, as shown in Fig. [8j peaks at lower values for stop 
quarks than for top quarks. An upper limit on Ht is 
effective at reducing the amount of tt background while 
a lower limit helps reduce the other backgrounds. 

A transverse mass cut is used in the low Am cut groups 
(a and b) to remove most of the DY background remain- 
ing after the other cuts have been applied. The transverse 
mass distribution for cut group b is shown in Fig. [9] 

The upper bound on the transverse momentum of the 
dilepton system is applied to the lowest Am cut group 
(a) to reduce the tt and diboson backgrounds. Its effec- 
tiveness decreases as Am increases. The sliding cut is 
determined by a linear fit to the optimal dilepton pt val- 
ues as a function of Am in cut group a: p 1 ^ < (cAm — 1) 
GeV/c. 



VII. RESULTS AND CONCLUSIONS 

The numbers of expected background and observed 
data events for cut groups a — d are given in Table PVTO In 
general the agreement between data and SM background 
estimations is good, although the statistical uncertain- 
ties are large. The largest deviation is an excess of data 
over background expectations for fifi events with large 
JSt- Examination of individual event properties found 
no evidence for cosmic rays or pion or kaon decays-in- 
flight. These results are not independent observations, 
since there is a large overlap in events between cut groups 
6, c, and d. The largest muon excess, found in cut group 



6, has a 0.4% probability for the modeled background to 
fluctuate up to the number of observed events or more 
(p- value), when considering this cut group alone. The 
p- value includes the effects of the estimated background 
uncertainty. Combining all channels in cut group b raises 
the probability to 12.4%. The corresponding one-sided 
Gaussian significances are 2.6 a and 1.2 a respectively. 

A joint likelihood is formed from the product of the 
individual channel likelihoods. Using this likelihood, we 
apply a Bayesian method [19] with a flat prior for the 
signal to set 95% confidence level upper limits on the 
production cross section at each considered point in the 
stop-sneutrino mass plane. Systematic uncertainties are 
incorporated by convolving the Poisson probabilities for 
the signal with Gaussian distributions representing each 
background uncertainty. The three dilepton flavor chan- 
nels are incorporated into the statistical analysis simul- 
taneously (but individually), with full treatment of cor- 
related and uncorrelated uncertainties. The correlations 
between the various backgrounds and the MC-generated 
signal are also accounted for. 

One-dimensional curves of the upper cross section lim- 
its and the theoretical cross section are shown in Fig. [10] 
for groups of points with fixed Am. The cross section up- 
per limits for a given Am tend to be almost independent 
of the stop mass. 

To determine the observed (expected) exclusion con- 
tour in the stop-sneutrino mass plane, we set the number 
of events equal to the number of events in the data (total 
expected background). We then calculate the 95% con- 
fidence level upper limits on the stop pair cross section 
for 74 points in the (m^, m^) plane to determine both 
the points we exclude and those that we expect to ex- 
clude in the absence of signal. We interpolate linearly 
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between nearby excluded and not-excluded points. The 
expected and observed limits are shown in Fig. [11] along 
with previous limits [H, ll| . 

In conclusion, this analysis extends the previously ex- 
isting exclusion limits to higher sneutrino masses for stop 
masses in the range 135-155 GeV/c 2 , and to stop masses 
up to 180 GeV/c 2 , for low sneutrino masses. For a par- 
ticular set of optimization cuts (cut group 6), a 2.6 a 
excess of fifi events is observed, but is reduced to 1.2 a 
when all channels are combined. 
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